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Abstract 

The purpose of this paper is to derive the anisotropic averaged Euler equa- 
tions and to study their geometric and analytic properties. These new equations 
involve the evolution of a mean velocity field and an advected symmetric tensor 
that captures the fluctuation effects. Besides the derivation of these equations, 
the new results in the paper are smoothness properties of the equations in 
material representation, which gives well-posedness of the equations, and the 
derivation of a corrector to the macroscopic velocity field. The numerical im- 
plementation and physical implications of this set of equations will be explored 
in other publications. 
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1 Introduction 



A fundamental problem in turbulent fluid dynamics is the difficulty in resolving the 
many spatial scales that are activated by the complicated nonlinear interactions. It 
is a challenge to produce models that capture the large scale flow, while correctly 
modeling the influence of the small scale dynamics. While there are many efforts 
in this direction, the goal of the present paper is to introduce a new method that 
is based on the combination of two basic ideas: the use of an ensemble averaging 
that represents a spatial sampling of material particles over small spatial scales, and 
the use of asymptotic expansions together with this averaging on the level of the 
variational principle. 

Our approach is conceptually similar to the method of Reynolds averaging and 
Large Eddy Simulation techniques, but has the advantage of 1) not needing an 
additional closure model and 2) automatically providing a small scale corrector to 
the macroscopic flow field. 

Our methodology has some interesting connections with the method of Optimal 



Prediction introduced by Chorin, Kast and Kupferman 199£], which will be explored 
in future publications. 

In the body of the paper we shall comment on a comparison between our ap- 
proach and that of Chen et al. 1998] and Holm I999|, which produces different 
equations. 



1.1 A Brief Review of the Euler and Isotropic Averaged Euler 
Equations 

A Brief History. There has been much recent interest in the averaged Euler 
equations for ideal fluid flow. In this paper we will focus on the geometry and 
analysis of a related set of equations, which we call the anisotropic averaged Euler 
equations. The original averaged Euler equations appear as a special isotropic case 
of the more general equations. 

The isotropic averaged Euler equations on all of M n first appeared in the context 



of an approximation to the Euler equations in Holm, Marsden and Ratiu [1998a | and 
some of its variational structure was developed in [Holm, Marsden and Ratm" 1998b] ; 
this variational structure retains the quadratic form of the variational structure for 
the original Euler equations, so that the equations can be viewed as describing a 



certain geodesic flow in a sense similar to the work of Arnold ]1966 ] and Ebin and 
T97Cfl . 

Remarkably, these equations are mathematically identical to the well-known 



Marsden 



inviscid second grade fluids equations introduced by Rivlin and Erickson 1955] ] . The 
geometric analysis of these equations, including well-posedness and other analytic 



properties, was developed in Shkoller [1998, 2000] and Marsden, Ratiu and Shkoller 
2000]. These references also discuss the relation to the second-grade fluid literature. 



In Oliver and Shkoller 200C |, the link with the vortex blob method was estab- 



lished; therein, it was shown that the vortex blob numerical algorithm generates 
unique global weak solutions to the averaged Euler equations. These weak solutions 
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induce a weak coadjoint action on the vector space of vorticity functions, modeled as 
the space of Radon measures. The existence of such a weak coadjoint action makes 
rigorous the formal constructions of Marsden and Weinstein 1983(| on the geometry 
of point-vortex and vortex blob dynamics. 



The works of Chen et al. 1998 1 and Holm 1999| formulated equations for the 



slow time dynamics of fluid motion by averaging over fast time fluctuations about 
the mean; that approach, founded on a Reynolds decomposition translated over 
the Lagrangian parcel, and the resulting system of equations, is different from our 
approach and from the results that we shall present. We give a few more details on 
the comparison in the body of the paper. 



The Euler Equations as Geodesies and Notation. It is well-known how to 
view the Euler equations as geodesies on the group of diffeomorphisms and that 



this view has concrete analytical advantages, due to the work of Arnold [1966 and 
Ebin and Marsden 1970(1 • I n particular, this work shows that the equations define a 
smooth vector field (a spray) on the group of diffeomorphisms, that is, in Lagrangian 
(or material) representation. The reduction of the equations from material to spatial 
(Eulerian) representation may be viewed by the classical and general technique of 
Euler-Poincare reduction (see Marsden and Ratiu 1999|| and Holm, Marsden and 



Ratiu 1998b] for an exposition and further references) and this view is a helpful 
guide to understanding other fluid theories as well. 

The geometric view of fluid mechanics, along with a careful understanding of 
the averaging process, will be basic to the present paper, so we briefly review the 
salient features of the theory for the reader's convenience, and to establish notation. 

Let (M,g) be a C°° compact, oriented n-dimensional Riemannian manifold with 
C°° boundary (possibly empty). Of course open regions with smooth boundary in 
the plane or space are key examples. The Riemannian volume form associated with 
the metric g is denoted fi. 

The Euler equations for the velocity field u of an ideal, incompressible, homoge- 
neous fluid moving on M (such as a region Q in IR 2 or M 3 ) are 



du 



+ (u ■ V)u = -Vp 



1.1] 



with the constraint div u = and the boundary condition that u is tangent to the 
boundary, dM. The pressure p is determined by the incompressibility constraint. 
The nonlinear term (u ■ V)u is interpreted in the context of manifolds to be V u u, 
the covariant derivative of u along u. In Euclidean coordinates, these equations are 
given as follows (using the summation convention for repeated indices): 



du 1 
~dt 



+ u J 



4 du 1 



dp 

dx l 



and on a Riemannian manifold (or in curvilinear coordinates in Euclidean space), 
the Euler equations take the following coordinate form: 
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where gij are the components of the Riemannian metric g, g %1 = [gij]~ , and T l - k 
are the associated Christoffel symbols. Using covariant derivative notation, these 
coordinate equations read 

We let the flow of the time dependent vector field u(t, x) be denoted by rj(t, x) so 
that 

d 

— r](t, x) = u(t,r](t,x)), 

with 7/(0, x) = x for all x in M. For each t, we denote the map rj(t, •) by r\t so that 
r/o = e, the identity map. Thus, the map x \— > r]t{x) gives the particle placement 
field for the fluid. Because of the incompressibility, each map r/t is volume preserving 
and is a diffeomorphism. 

We shall be working with vector fields u of Sobolev class H s for s > (n/2) + 1 
and, correspondingly, r/t £ T>^, the group of fP-volume preserving diffeomorphisms. 
If there is any danger of confusion, we shall write T> ^(M) to indica te the underly- 
ing manifold M. See Ebin and Marsden 1970 1 and j5hkollcr~ 2000| for some basic 



properties of Hilbert class diffeomorphism groups for manifolds with boundary. 

Arnold's theorem on the Euler equations may be stated as follows: A time 
dependent velocity field u satisfies the Euler equations iff the curve r\t is a geodesic 
of the right invariant L 2 -metric on . 

This L 2 -metric is defined as follows. The tangent space to T>^ at the identity is 
identified with the space X^ iv , the space of H s divergence free vector fields on M 
that are tangent to the boundary dM. The right invariant L 2 -metric is defined to 
be the weak Riemannian metric on T> s whose value at the identity is 



(u,w) L 2= {u(x),w(x)) xfJ ,(x), 
J M 

where we write the pointwise inner product as (u(x),w(x)} x = g(x)(u(x),w(x)), 
and the pointwise norm |tt(2;)| 2 = (u(x),u(x)) x . 

As we shall explain shortly, with the maturation of Euler- Poincare theory, Arnold's 
theorem becomes an easy consequence of more general and rather simple results. 

Lie Derivative and Vorticity Form. As is well-known, the Euler equations can 
be written in terms of Lie derivatives as 

^ + £^ = dQ|n| 2 -p) =-dp>, (1.2) 

where u is the one-form associated to the vector field u via the metric, and 
denotes the Lie derivative of the one-form u along u. Taking the exterior derivative 
of (|L^) gives the familiar advection equation for vorticity: 

duo 

— + £ u uj = 0, 
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where uj = du^ is the vorticity, thought of as a two- form. In 2D, uj is identified with 
a scalar and is traditionally thought of as the 2D-curl of the velocity field, while 
in 3D, uj may be identified (using the volume-form jj) with a vector field which is 
traditionally obtained by taking the curl of u. 

The vorticity equation is the infinitesimal version of the following advection 
property: 

&t = {Vt)*vo- 

Of course in two dimensions, this gives the usual advection of vorticity as a func- 
tion, while in three (or higher) dimensions, the advection is understood in terms of 
advection of two- forms. 

The Euler equations have both an interesting Hamiltonian structure in terms 
of Poisson brackets (a Lie-Poisson bracket) and a variational structure. In this 
paper we shall be working primarily with the variational structure; the Hamiltonian 



structure, along with references to the literature may be found in |Marsden and 



Weinstein [1983| , [Arnold and Khesin [1998| and [Marsden and Ratiu [1999 ] 



Lagrangian and Variational Form. The Lagrangian is given by the total kinetic 
energy of the fluid; in spatial representation, this Lagrangian is 

L{u) = \ [ |n(x)|V (1-3) 

1 JM 

The corresponding (unreduced) Lagrangian on TP* is given by 



£(v,v) = t; [ g{v(x)){v(x),v(x))v- (1-4) 

' M 



2 



Hamilton's principle on D* applied to the Lagrangian C gives geodesies on this 
group. Euler-Poincare reduction techniques (see Marsden and Ratiu [1999| ) show 



that this variational principle reduces to the following principle in terms of Eulerian 
velocities: 

5 I L(u) dt = 0, 

J a 

which should hold for all variations 5u of the form 

Su = w + [u, w], 

where to is a time dependent vector field (representing the infinitesimal particle dis- 
placement) vanishing at the temporal endpointsQ. Here, [w,u] denotes the usual 
Jacobi-Lie bracket of vector fields. One readily checks that this reduced princi- 
ple yields the standard Euler equations. This simple computation is the heart of 
Arnold's theorem. 



The constraints on the allowed variations of the fluid velocity field are commonly known as "Lin 
constraints". This itself has an interesting history, going back to Ehrenfest, Boltzmann, Clebsch, 
Newcomb and Bretherton, but there was little if any contact with the heritage of Lie and Poincarc 
on the subject. 
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Analytical Issues. While the Eulerian (spatial) representation has been empha- 
sized in most analytic studies of the Euler equations, fluid motion viewed on the 
Lagrangian (material) side has some distinct advantages. For example, [Ebin and 
Marsden [1970 1 proved that the flow, solving the Euler equations, on the volume- 
preserving diffeomorphism group T>^, s > n/2 + 1, is smooth in time. They derived 
a number of interesting consequences from this result, including theorems on the 
convergence of solutions of the Navier-Stokes equations to solutions of the Euler 



equations as the viscosity goes to zero when M has no boundary. In addition, Mar- 



chioro and Pulvirenti [1994 1 analyzed the Lagrangian flow map to establish sharp 
well-posedness of the 2D Euler equations and prove convergence of the vortex blob 
algorithm. In many cases, the Lagrangian framework is, in fact, the more natural 
setting to study the behavior of solutions, and we shall emphasize this point of view. 



1.2 The Averaged Euler Equations 



The Isotropic Averaged Euler Equations. Let a be a positive constant. In 
Euclidean space and in Euclidean coordinates, the isotropic averaged Euler equations 
(inviscid second-grade fluids equations )^ read: 



dt dxi 



a 



An,- 



dp 

dx l 



where v = u — a 2 Au and A denotes the componentwise Laplacian, and there is a 
summation over repeated indices (in Euclidean coordinates, as is common, we make 
no distinction between indices up or down). While there are several choices, the no 
slip boundary conditions u = are often used for this model. 



Rate of Deformation Tensor. One of the interesting things that comes out 
of a careful derivation of the equations is the natural occurence of the rate of 
deformation tensor, which is defined by 

Defu = - (Vu + (Vu) T ) 
which we write in coordinates as: 

(i>«r»)' ; ^(» ; : , • «',)■ 

We also let Def -u b = \ [Vu b + (Vu b ) T ] which we write in coordinates as 

D i:j = ( Def v 9 ) .. = - [m-j + u j;i ) . 

Note that this is exactly the Lie derivative of the metric tensor; that is, Def u b = £ u g, 
which is sometimes called the Killing tensor. 

2 These are also known as the Euler-a equations. 
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Smoothness Properties. Results on smoothness of the Lagrangian flow map for 



given in 


Shkoller 


199^1 on compact boundary- 


Marsden, Ratiu and Shkoller 


2000 on compact 



Euclidean domains. The problem of how to formulate this system on compact Ric- 



mannian manifolds with boundary was solved in Shkoller 200C | ; the equations take 
the form 



d t (l - a 2 A r )u + V u (l - a 2 A r ) u — a 2 (Vii)* • A r u 



gradp, 



together with the constraint div u = 0, and with appropriate initial conditions u(0) = 
uq, as well as boundary conditions. The symbol A r is the operator Def* Def acting 
on divergence- free vector fields, where Def* is the L 2 formal adjoint of the (rate of) 
deformation operator Def. Explicitly, 



A,. 



-(dS + 5d) + 2Ric. 



;i.5) 



As with the usual Euler equations, the function p is determined from the incom- 
pressibility condition. 

Lie Derivative Form — The Isotropic Equations. The averaged Euler equa- 
tion can be neatly written in terms of Lie derivatives: 

d t v b + £j 



-dp, 



(1.6) 



where v" = (1 — a 2 A r )u b . 

The Anisotropic Averaged Euler Equations. These equations, which are the 



main subject of the present paper, and which will be derived in §^2, will now be 
stated. The basic variables that are evolving in the anisotropic averaged Euler 
equations are the macroscopic velocity field u and a symmetric tensor field F on 
M; the tensor field F will be interpreted as the contravariant spatial fluctuation 
tensor and it will keep track of the anisotropy of the fluid deviations from the 
macroscopic flow. These equations also depend on a choice of length scale a. 
It is convenient to define the linear operator C : X^ iv n Hq — > H s ~ 2 , s > 1, by 

Cu := Div \c : Vu b ] , 

where b is the map from vector fields to one-forms associated with the metric g, and 
the fourth-rank symmetric positive tensor C is the symmetrization of the tensor 
F (g) g~ l , given in local coordinates by 



C 13 



rikl 



I (W fe + f kj 9 U + F li g jk + F fc V ') • 



With this notation, the anisotropic averaged Euler equations on manifolds 



arc 



d t (l - a 2 C)u + V u (l - a 2 C)u - Q 2 [Vuf • Cu + 2a 2 F : V(Def rf) 



>\2 



4a 2 Div ( (Def ti) 2 • F 



gradp, 
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together with the advection equation 

d t F + £ U F = 0, 

the incompressibility constraint div u = 0, initial data u(0) = uo and F(0) = Fq, 
and the Dirichlet boundary condition u = 0. 



Lie Derivative Form — Anisotropic Equations. The anisotropic averaged Eu- 
ler equations can also be written using Lie derivatives as 



d t v b + £ u v h + 



2a 2 F : 



V(Def v°) 2 - Aa z Div ( (Def uf ■ F 



-dp, (1.7) 



where v" = (1 - a 2 C b )u b , where CV = Div[C : Viff. 



Coordinate Form. In local coordinates, the anisotropic averaged Euler equations 
become 

dt (u* - a 2 [C^ kl n k ,j U ) + ("*- A&^k* },l ) m u m - a 2 u m)l [C m ^ l u k ,j U 
= 

together with the advection equation 

d t F^ + F^, k u k - F k ^u\ k -F ik u\ k = 0, 

with the constraint u l i = 0, given initial conditions u l (0) = u l , and with the no-slip 
conditions u l = on the boundary. If the metric g is not the Euclidean metric 
5ij, then the partial derivatives above should be interpreted as arising from the 
Levi-Civita covariant derivative associated to g. 



1.3 Outline of the Main Results. 

The main results of the present work are as follows: 

1. We derive, in a systematic way, the first order averaged Lagrangian given in 
coordinates by 

Lf(u, F ) = IJ m {9iku l u k + 2a 2 g ik Fi l D ij D kl } [det g^dx. 

and, using the calculus of variations, derive the associated anisotropic averaged 
Euler equations as the corresponding Euler-Poincare equations. The Euler- 
Poincare technique was also used in Holm 1999(| , but the Lagrangian and 
associated equations are different. In particular, the principles and philosophy 
governing the derivation of the Lagrangian are completely different. 
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We show that the equations are well posed; in fact, we show more, namely that 
the corresponding Lagrangian flow map is smooth in time in the appropriate 
Sobolev topology. 

Another important achievement is that while the macroscopic velocity field u 
is computed on spatial scales larger than a, we are able to obtain a corrector 
for this macroscopic field to order a 2 . This is done in §4.2 and is similar to 
what one does in the theory of homogenization. 



2 The Derivation 
2.1 Introduction 

This section presents a new method for constructing models of hydrodynamics which 
takes into account the fundamental idea that a fluid particle is not a point, but rather 
a collection of points forming a representative sampling. Our approach is founded 
upon a certain type of Lagrangian ensemble averaging performed at the level of the 
variational principle. A similar idea on the level of the equation itself, as opposed 
to the variational principle, was used by Barenblatt and Prostokishin 1993| for 
deriving models of damage propagation. 



Naive Averaging Does not Work. We first explain why the naive approach to 
spatially averaging a quadratic Lagrangian or Hamiltonian does not suffice. As a 
simple example, consider the Lagrangian on scalar functions on M n given by L(u) = 
\ Iw n u ( x )dx and for a given positive constant a, define a new averaged Lagrangian 
by 

L a (u) = [ / / u\z)dzdx 

Jm« \rS{X,a)\ JB{x,a) 

which is obtained from L by averaging the original Lagrangian over balls of radius 
a. Here B(x, a) denotes the ball of radius a about the point x in M n and \B(x, a)\ 
denotes its volume. 

Taylor expanding the integrand about x and then integrating by parts yields 
cancellation of all but the zeroth-order term, thus reproducing exactly the original 
Lagrangian L. This is to be expected since the quadratic nonlinearity is rather 
weak, and since absolutely no information concerning the local spatial structure of 
the continuum is being provided. The latter issue is of fundamental importance and 
is the foundation upon which we shall build our theory. 



2.2 The Averaging Construction. 

To implement our construction, we will average over an ensemble of Lagrangian 
fluctuation maps. We will now proceed to develop this formalism. 



2.2 The Averaging Construction. 10 



Fuzzying the Lagrangian Flow. Let (M, g) be a C°° compact, oriented re- 
dimensional Riemannian manifold with C°° boundary (possibly empty). We con- 
sider a two-parameter family of volume-preserving diffeomorphisms of M de- 
pending on a "radial" component e G [—a/2, a/2], a > 0, and an "angular" compo- 
nent 9 € <S"~ , where denotes the upper hemisphere of the unit sphere S 1 ™" 1 
in W 1 . In case M has nonempty boundary, we embed M into its double M and 
consider this two-parameter family defined on M; in this case, £ e ' need not leave 
dM invariant. This fact will be important later for certain ellipticity properties. 
The parameterization is chosen such that 

£°> e {t,x) = x, 
dist(x,£ € - e (:z;)) < |e| 

for all e S [—a/2, a/2], all t, and 9 € We define the infinitesimal fluctua- 

tion vector by 

a vector field depending on the parameter and time t. 

For each time t, the Lagrangian flow map rjt, where rft(x) = rj(t,x), associated 
with a solution of the Euler equations is a volume-preserving diffeomorphism of M 
which maps fluid particles x € M to rjt(x) € M. Motivated by the idea that a 
particle in a continuum should really be regarded as a representative of a sample of 
particles over a region, we define the ^'^-perturbed particle placement field by 

4'\x) = (Ct 6 r 1 °r }t {x) (2.1) 

for all e € [—a/2, a/2] and G 5"™ _1 . The family of maps n|' is called the fuzzy 
flow. For each e, and £, the map : M — > M is a volume-preserving diffeomor- 
phism of the fluid container. Note that at e = 0, t^' 9 = rjt for all 9 € S 1 " -1 . 



We take ^ e X>*(M) and ^' G ^ ( M ) so that r?^ G D"(Af). See §L1| for the 



definition of the group P 



Decomposition of the Spatial Velocity Field. Our goal is to derive the Eu- 
lerian velocity field uf corresponding to the £^' 9 -perturbed particle placement field 
rf t ,d , and define a new Lagrangian by averaging the velocity u t ,e over the radial 
parameter e and the angular coordinate 9. We shall proceed with this averaging 
process as follows: we begin by defining the Eulerian vector fields associated with 
our three Lagrangian maps. Let 

d t r](t,x) = u(t,rj(t,x)), 

d t C' e (t,x)= W ^(t,e fi (t,x)), 

d t 7] e ' e (t,x)=U^(t, V e ' 6 (t,x)). 

Differentiating the Lagrangian decomposition ( |2.1| ) with respect to time t, we 
obtain the spatial velocity decomposition 

u^(t,x) = {e> e *(u-w^)}(t,x), (2.2) 
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where the notation £ £ > e * denotes the pullback by the map £ e ' e . We can also write this 
decomposition using the push-forward notation via the relation (d' 6 )* = (Ci'^) 
so that the action on a vector field v is given by 



(e/yv = T(e t ' e r i ovoc i 



where we use the symbol T to denote the tangent map (which is locally represented 
by the matrix of partial derivatives). Thus, the decomposition (|2.2j ) may be equiv- 
alently written as 

u^(t, X ) = ne/rHx) o (u( t ,e t > e (x)) - «,«■«(*, #«(*))) , (2.3) 

where, again, u e ' (t,x) is the Eulerian spatial velocity field corresponding to the 
fuzzy flow r] t ' . 

Comments on the Nature of the Decomposition. The Lagrangian decompo- 
sition (^^) which "fuzzies" the Lagrangian flow map yields the decomposition fl2.3| ) 
for the corresponding Eulerian variables which is of a hybrid Lagrangian- Eulerian 
type. The Lagrangian characteristics of this decomposition are encompassed in the 
presence of the purely Lagrangian fluctuation maps £ t ' , and it is indeed the pres- 
ence of this Lagrangian term in (^^) which allows us to proceed with an asymp- 
totic expansion which is both philosophically and mathematically different from the 
"naive" expansion we discussed earlier. We should also emphasize that without 
this Lagrangian aspect, the decomposition ( [2.3| ) would reduce to the usual additive 
(Reynolds) decomposition of spatial velocity fields into their mean and fluctuating 
parts, which does not reflect the fuzzyness of the Lagrangian flow map. 

Our approach should also be contrasted with the approach taken by phen et~aT 
19981 an d |Holm |1999| . In those papers, the decomposition 



rf{t, x) = n(t, x) + er(f, rj(t, x)), 

is made, where o~(t, x) is a fluctuation vector field, and rj a (t,x) is a perturbed 
Lagrangian trajectory of the reference element x. This decomposition is intrinsi- 
cally problematic, in that a material vector field a(t,r)(t,-)) is being added to a 
volume-preserving diffeomorphism r](t, •). As a consequence, the perturbed trajec- 
tory r] a (t,x) does not come from a volume-preserving diffeomorphism of the fluid 
container, that is, rf(t, •) is not a volume-preserving map. 

The Averaged Lagrangian. We define the averaged Lagrangian L a by 

L a (u) = \ [ - f (u*> e (t,x),u*> e (t,x))de v(6) fi(x) 

2 J M « J -a/2 J si- 1 



1 r 1 ri 2 

2 JM a J -a/2 Js$- 



e> 6 *(u-w e > e )(t,x) 2 ale i/(0) fi(x), 



3 This decomposition can also be written as if' (t,x) = Ad, t e.e, x (u t — wl' ), where Ad is the 
adjoint action of the volume-preserving diffeomorphism group on divergence-free vector fields. 
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where /j, is the Riemannian volume form on M, and v is the induced Riemannian 
volume form on the upper hemisphere of the unit sphere S n ~ x in R n . 



Comments on the Nature of the Fuzzying Operation. By using the upper 
hemisphere S"™ -1 and integrating from —a/2 to a/2, we are tacitly assuming that 
there is a hyperplane of symmetry in the ^-parameter space. This is not a restric- 
tion even near the boundary, as the hyperplane of symmetry can always be chosen 
orthogonal to the boundary and the maps £ t ' can be chosen to be symmetric about 
this hyperplane with respect to the radial parameter e. 

The reader should keep in mind that the variables 8 and e parameterize possible 
families of maps and are not to be confused with spatial spheres in the flow itself. 
We are averaging over these families of maps and not literally over spatial regions. 
A representative of the family of fluctuation maps in the two dimensional case 



and near a boundary is shown in Figure 2.1 





Figure 2.1: An example of a perturbing map its inverse takes the flow point rjt(x) to the 



perturbed flow point rf t ,e (x). The parameter space for (e, 6) and the symmetry plane is shown on 
the left. 

The internal structure behind the fuzzyness of the macroscopic Lagrangian 
is completely encoded in the fluctuation maps £ t ' . 

The zeroth-order assumption that these maps are simply the identity map leads 
to the zeroth-order Lagrangian Lq which is exactly equal to the Lagrangian L 



given in (1.3) and thus produces the usual Euler equations of hydrodynamics as the 
continuum model. We proceed to obtain the first order correction to this model 
which accounts for the spatial fluctuations. 



4 For a general continuum, the information about the structure of the representative sampling 
would be encoded in the fluctuation maps. For example, this might include defects or microstruc- 
ture. 
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Asymptotic Expansion. We Taylor expand u e ' in e about e = to obtain 

u e ' e (t, x) = Ct' e *{ u t ~ w e t ' 9 )(x) = u t (x) + €£^ e>t)X )Ut(x) - ei'(6, t, x) + 0(e 2 ), 

(2.4) 

where the overdot means the time derivative. This follows from the definition of 
the Lie derivative, the fact that w e,s = dt^ ,d o , and that ^ = e. Using the 

zero-torsion condition on the Levi-Civita covariant derivative, £^u = V £/u — V u £', 
and suppressing the dependence on t and x, we get 

u € ' e = u + e(V« • £{6) - Vf'(0) • u - i'{6)) + 0(e 2 ) 

or, in index notation, 

u ei = u l + e (u l .fi{6) - i' :j {0)u-< - + 0(e 2 ), 

where 

u = u i d i , ? = C9i and Vu • f = u'.jg A. 

In order to proceed, we make the first-order Taylor Hypothesis that the 
infinitesimal fluctuation vector £' is frozen, as a one- form, into the fluid so that its 
Lie transport vanishes; namely, 

(£) b + A(O b = 0. (2.5) 

We again express the Lie derivative of the 1-form field (£') in terms of the covariant 
derivative to obtain, in index notation, 

ii + ufa + = 0. 

From this hypothesis, the O(e) term in the Taylor expansion Q2.4] ) is 

= Ui- d e + Uj -4 j 2l)of u'-i'(O). 

It follows that 



u 



efi __ ^,6y(u- w *)= u + 2e'Defu-£'(e) + 0(e 2 ). (2.6) 



L a {u) = \ I - I I" 2 [Kx)| 2 + 2e(n(x),Defn(x) • £'(s,0)) 
2 JM « Js"- 1 J -a/2 



Substitution of (2.6) into the averaged Lagrangian L a yields 

'■a/2 

2 Jm a Js^ 1 J -a/2 

+ 4e 2 |Defu(x) ■ £'{x,9)\ 2 + 0{e 3 )} de v(6) fi(x). (2.7) 

An important point about this expression is the following: There is no contribu- 
tion from the term (u,0(e 2 )} to the energy at order 0(e 2 ). In fact, the 0(e 2 ) term 

in has the form 0(e 2 ) = e 2 ((" + R) , where 
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However, £" is an independent field and must have its own dynamics specified. We 
assume that this dynamics is chosen so that £" + R is O(e) and so the a priori 
0(e 2 )-term in fl2.6|) is in fact 0(e 3 )- 

Integrating ( |2.7D in e, rescaling a i— > Wa/6, and defining the symmetric rank-2 
contravariant spatial fluctuation tensor (indices up) F by 



f{x) = / { , (i 1 e)«f'(s,eM9) I 

is™- 1 

we obtain the first- order averaged Lagrangian 

L?(u,F) = l[ [ [\u\ 2 + 2a 2 \Defu-£'{x,9)\ 2 ] u{9) a{x) 

= - [ [\u(x)\ 2 + 2a 2 {F(x) o Def ^(x), Def u(x))] u(x) 
2 Jm 

= o/ {l«(^)| 2 + 2a 2 [g(x) ® F(x)] : [Def u(s) ® Def u(x)]}^(x). (2.8) 
2 Jm 

In coordinate notation, the first-order averaged Lagrangian takes the form 

L?(u, F ) = \f M {^ku l u k + 2a 2 g lk Fi l [Def u^Def it],* } [det ^adx. 

The first-order averaged Lagrangian is a function of the macroscopic Eulerian 
velocity field u and the contravariant spatial fluctuation tensor F. 

The Isotropic Case. In the case that the fluctuation tensor is isotropic so that 

F(x) =g~\x), 

the isotropic first- order averaged Lagrangian L^ iso is given by 



LlUu) = \ [ [\u(x)\ 2 + 2a 2 \Defu(x)\ 2 ]u(x). 
1 ->M 



In this special case, the Lagrangian depends only on the Eulerian velocity field u 
and no semi-direct product theory is required; in fact, the standard Euler-Poincare 
theory for reduced Lagrangian variational principals may be invoked to obtain the 
isotropic averaged Euler equations as 

<9 t (l - a 2 A r )u + V u (l - a 2 A r )u - a 2 (Vii)' • A r u = - gradp, 
div u = 0, u(0) = uo, 



where A r = —{d5 + Sd) + 2Ric (see [Bhkoller [2000| ) . As we stated above, the 



equations (|2.9| ) are precisely the equations of inviscid second-grade non-Newtonian 
fluids, and exactly coincide with Chorin's vortex blob algorithm when a particular 
choice of smoothing kernel is used (see Oliver and Shkoller 2000|j ). In the case 



that M is a manifold without boundary, the incompressibility of the fluid allows us 
to replace the term 2a 2 1 Def u\ 2 with a 2 |V«| 2 in (|2.8| ), and still obtain the identical 
evolution equations as in (^1]); however, for domains with boundary it is essential to 
retain the strain tensor Def u in the Lagrangian so as to obtain the natural boundary 
conditions which ensure ellipticity of the operator (1 — a 2 A r ). 



3 The Variational Principle and Semidirect products 
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3 The Variational Principle and Semidirect products 

We shall next explain the sense in which the Lagrangian Lf(u,F) defined in ( p.8[) , 
a function of spatial variables, can be obtained from a Lagrangian defined in ma- 
terial variables. This will be done via an Euler-Poincare procedure, which involves 
the group T> s D (§) C°° (T 2,0 (M)) , the semi-direct product of the volume-preserving 
difieomorphism group T>^ D (with Dirichlet boundary conditions) and the smooth 
sections of the vector bundle T 2 '°(M), consisting of second-rank contravariant sym- 
metric tensors. Before proceeding to our specific example, we shall digress briefly 
to explain the general theory. 

3.1 Lagrangian Semidirect Product Theory. 

The General Set Up. Let V be a vector space and assume that the Lie group 
G acts linearly on the right on V (and hence G also acts on its dual space V*). In 
the case that the vector space V consists of sections of a vector bundle E, V* will 
denote the sections of the dual bundle E*. The semidirect product S = G®V is 
the Cartesian product S = G x V whose group multiplication is given by 

(gi,vi){g2,v 2 ) = (gig2,v 2 + vig 2 ), (3.1) 

where the action of g G G on v € V is denoted simply as vg. The Lie algebra of S 
is the semidirect product Lie algebra, s = g(S)y, whose bracket is 

[(V>i,«i), (ip 2 , v 2 )] = ([ipx,tp 2 ], viil) 2 - v 2 ipi), (3.2) 

where we denote the induced action of g on V by concatenation, as in v 2 ipi- For 
v € V and a £ V* , define the bilinear operator v o a € g* by 

(v o a , 77) = (arj, v) . 

The Objects in Our Case. We choose G to be the topological group *D S D . 
While this is not a Lie group, right multiplication is a smooth operation, and this is 
the crucial feature we shall make use of. The tangent space at the identity T e T> s ^ D 
is equal to 3tj iv D , the H s vector fields on M vanishing on the boundary and with 
zero divergence, and plays the role of the Lie algebra g. 

We set V = H S (T 2 '°(M)), the H s sections of the vector bundle T 2 '°(M) consist- 
ing of contravariant symmetric two tensors (indices down). Thus, the vector space 
V* is H S (T 0,2 (M)), the H s sections of the covariant two tensors (indices up). The 
duality is with respect to the induced Riemannian metric on T 2,0 (M). The topo- 
logical group T> s D acts on the vector space H S (T 2,0 (M)) by pull-back; hence, this 
action takes values in H B ~ 1 (T 2 >°(M)). Since the group is volume preserving, the 
induced right action on V* is also by pull-back. We have the map (r],F) 1— > r]*F. 

It follows that the infinitesimal action is by the Lie derivative which also maps 
H s sections into sections of class H s ~ 1 . Thus, according to the above definition, 



3.1 Lagrangian Semidirect Product Theory. 
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the diamond operator is computed as follows: Let K 6 H s 1 (T 2 '°(M)), F G V* = 
H S (T°' 2 (M)) and let u G Q = X s div D . We define the operator 

£f ■ n*,D - H s ~ 1 (T 2 '°(M)), = 

Then the adjoint operator (with respect to the Riemannian metric on H S (T 2,0 (M))) 
C* F : H s ~ l {T^ 2 {M)) -» 3^ ivD and is defined by 

(A' oF,u) = (£ U F, K) = (u, C* F K) . 

Thus, we have 

K o F = C* F K. 

Semidirect Euler-Poincare Reduction. Assume we have a right G-invariant 
function £ : TG x V* -» R. For a G V*, let £ ao : TG -> R be given by C aa (v g ) = 
C(v g ,ao), so £ ao is right invariant under the lift to TG of the right action of G ao on 
G, where G ao is the isotropy group of ao- Define L : q x V* — ► R by 

L(ug5 _1 ,ao5 _1 ) = C(v g ,a ). 

For a curve <?(t) G G, let £(t) := g(t)<jf(t) _1 and let a(t) = aog(i)" 1 , which is the 
unique solution of the equation a(t) = —a(t)£(t) with initial condition a(0) = ao- 
In our setting, a = F G G°°(T 2 ' (M)), 

^ = u„ G {v G # S (M, TM) | v o G £^ iv , D D H^(TM), V G 2^ D }, 

and L(u v o r]~ 1 ,i]*Fo) = C(u v ,Fq) where L is given by Q2.8| ). 

Theorem 3.1. TTie following are equivalent: 

i Hamilton's variational principle 

5 I " C ao (g(t),g(t))dt = (3.3) 

holds, for variations b~g{t) of g(t) vanishing at the endpoints. 
ii g(t) satisfies the Euler-Lagrange equations for C aa on G. 
iii The constrained variational principle 

5 [ 2 L(£(t),a(t))dt = (3.4) 

holds on q x V* , using variations of the form 

ti£ = V-[€,v]i Sa = -an, (3.5) 
where rj(t) G vanishes at the endpoints. 
iv The Euler-Poincare equations hold on g x V* 
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3.2 Computation of the Anisotropic Averaged Euler Equations 

It is convenient to define the linear operator C : X^ iv Hifp — > H s ~ 2 , s > 1, mapping 
divergence- free vector fields to vector fields, by 



Cu := Div 



C : Vu b 



where b is the map from vector fields to one- forms associated with the metric g, 
and again the fourth-rank symmetric positive tensor C is the symmetrization of the 
tensor F ® g~ l , given in local coordinates by 

The functional derivatives of with respect to it and F are given by 

6Ll - (1 - a 2 C)u 



and 

We can then compute that 

6L " oF = 2a 2 F 



8u 



^ = 2a 2 [Defu] 2 . 



V(Def «' 



\>\2 



SF 

Letting t = (Def u) 2 , in index notation, we get 
'5Lf 



4a 2 Div ((Def u) 2 ■ F 



5F 



oF 



2o 2 /••'•'/,,;/, - 4a 2 [F ij t 



Using Theorem 3.1, we derive the following result. 

Theorem 3.2. The Euler- Poincare equations on Riemannian manifolds, associated 
to the Lagrangian L" given by (ft-Sj) , are the following anisotropic averaged Euler 
equations: 



h2 



d t (l - a 2 C)u + V u (l - orC)u - a 2 [Vuf ■ Cu + 2a 2 F : V(Def it") 



- 4a 2 Div ((Def u) ■ Fj = - gradp 
together with the advection equation 

8 t F + £ U F = 0, 



(3.7) 



(3.8) 



the incompressibility constraint div u = 0, initial data u(0) = uq and F(0) = Fq, 
and no-slip conditions u = on the boundary. 
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Anisotropic Averaged Euler Equations in General Coordinates. In gen- 
eral coordinates on a manifold, the averaged Euler equations read 

d t (u* - a 2 [C^ kl u k ,j],i ) + («*- AC ijU u k ,iU ) , m u m - a 2 Umn [C m i kl Uk ,j},i 
+ 2a 2 F^[D km g mn D nJ U -4a 2 [F k W im g mn D nj ], k 
= ~P,i 

where, as earlier, D{j = i (m-j + Uj-j) is the rate of deformation tensor and indices 
are raised and lowered using the metric tensor (which of course need not be diag- 
onal in general coordinates), and C ijkl = \ (F lj g lk + f kj g d + F h g jk + F kl g jl ). In 
Euclidean space, one need only set the components of the metric tensor gij to the 
Kronecker delta 

Comments on the Form of the Equations. In 2D, identifying F with the 
vector (F u , F 12 , F 22 ), equation (Kjfy takes the form 



D 



- F n - 




F 12 




p 22 





2u\ 



2u) 2 


6 2u 2 







u 



,2 

-2U 1 





- F n - 




pl2 




F 22 



(3.9) 



where D /dt denotes dt + (u ■ V). Notice that the matrix on the right-hand-side of 
(|3.9| ) is traceless; a similar form holds in 3D as well. This is not surprising, since by 
virtue of the incompressibility of the Lagrangian flow and the fact that Ft = f?t -^b, 
we have that 

det(F 4 ) = det(Fo), 

for all t for which the solution exists. As consequence, the operator (1— a 2 C) remains 
uniformly elliptic, if Fq is strictly positive. 



The Circulation Theorem. Let 7 : S 1 — ► M be a loop and let 7* = r/t 07 denote 
the evolution of the loop moving with the fluid. 

Theorem 3.3. For a solution of the anisotropic averaged Euler equations, we have 



2 r \> 

a C 



2a z 



it 



2 Div ( (Def u) ■ F) — F : V(Defti 



This follows directly from the Lie derivative form of the equations given in ( |Li 
We note that if one were to make use of the general Kelvin-Noether theorem given 
Holm, Marsden and Ratiu 1998b], one would arrive at the same result. 
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4 Analytic Properties 

In this section we prove well-posedness and other properties of the solutions by 
showing that these equations are given by a smooth vector field in material repre- 
sentation in the appropriate Sobolev topologies. This is in line with what is known 
about the Euler equations, as described in the introduction. We also discuss the 
corrector for the equations. 



4.1 Well-posedness of Classical Solutions 
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4.1 Well-posedness of Classical Solutions 

We shall prove existence, uniqueness, and smooth dependence on initial data on 
finite time intervals for solutions of the anisotropic averaged Euler equations. For 
simplicity, we shall restrict the fluid domain M to be a compact subset of Euclidean 
space with smooth boundary, although our methods can be applied to Riemannian 
manifolds. 

We begin by collecting some preliminary results. Set V s = Hq n H s and = 
Hq nX^ iv . Also, let T> S D denote the H s class diffeomorphisms which fix the boundary, 
and again let D denote the diffeomorphisms in T> S D which preserve the volume \i. 

Lemma 4.1. For u E V?,, s > 1, 

d t Cu = C{d t u) + Div [-V U F ■ Vu + Vu • Vu • F + Vu • F ■ Vu 1 ) , 

and 

V u Cu = C(V u u) + Vu • V u Div F + VVu : V U F - 2 VVu : (Vu • F) 
- Vu ■ (Vu • DivF) - Vu • (VVu : F). 

Proof. The proof is a simple computation which we leave to the interested reader, 
c.f. Lemma 3 in [Shkoller [200C|l . ■ 

Set C = Def* [(g ® F) : Dei]. Then C is a positive unbounded self-adjoint oper- 
ator on L 2 with domain V 2 . Define the inner-product (•, •) on V 2 by 



(u,v) = ((1 — a 2 C)u,v) 



L- • 



For s > n/2 + 1, (•,•) defines an inner-product on T e X>* the tangent space at the 
identity of the subgroup T> s ^ D consisting of those elements of which restrict to 
the identity on dM. Right-translating (■,■) to the entire group D defines a C°° 
weak Riemannian metric by Proposition 3 of Shkoller 2000| . 

Proposition 4.2. For r > 1 we have the following well-defined decomposition 

V = Vl l ®(l-£)~ 1 gradH r - 1 (M). (4.1) 

Thus, if F G V r , then there exists (v,p) eVJx iT r_1 (M)/R such that 

F = v + (1 — gradp 

and the pair (v,p) are solutions of the Stokes problem 

(1 - £)u + gradp = (1 - C)F, 

divu = 0, (4.2) 
v = on dM. 

The summands in are (•, •)- orthogonal. Now, define the Stokes projector 

P e (F) = F — (1 — Z:)" 1 grad p. 1 j 



4.1 Well-posedness of Classical Solutions 
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Then, for s > (n/2) + I, V : TT> S D — > TX>* D , giwen on eac/i ./i&er fry 



V v (X r] )=[P e (X v or ] - 1 )] or,, 
is a C°° bundle map covering the identity. 

Proof. The proof is identical to the proof of Proposition 2 in Shkoller [2000(1 • 



Theorem 4.3. Set s > (n/2) + 2, and let ((•, •)) denote the right invariant metric 
on T>^ D given at the identity by (•, •). For uq € T e T>^ D and Fq € C^iT 2 ' ), there 
exists an interval I = (— T, T), depending on \uq\ s , and a unique geodesic r] of ((•, •)) 
with initial data i](0) = e and ?)(0) = no such that 

Vis in C°°(I,TV^ D ) 

and has C°° dependence on the initial velocity uq. 

The geodesic rj is the Lagrangian flow of the time- dependent vector field u{t, x) 
given by 

d t n(t,x) = u(t,r](t,x)), 
and, with F(t,x) = (rj t ) # Fq(x) , 

(u,F) e C°(I,VP HCHiy^ 1 ) x C°(I,H S -\T 2 ' )). 

uniquely solves the anisotropic averaged Euler equations with Dirichlet boundary 
conditions u = 0, and depends continuously on (uq,Fq). 

Proof. The key to the proof rests in the fact that the pair (n, F) solves the anisotropic 
averaged Euler equations if and only if n is a solution of 

fi + U a (r) o if 1 ) o i] = [(1 - a 2 Cy l gradp] o n, (4.4) 

where 

U a {u) =a 2 (l - a 2 LY x {Div [V U F • Vn — Vn • Vn • F — Vn • F ■ Vn*] 
- Vn • (V« Div F) - VVn : V U F + 2VVn : (Vn • F) 
+ Vn • (Vn • Div F) + Vn • (VVn : F) - Vn* • Cu 
+2F : V(Def n 2 ) - 4 Div [F • Def n 2 ]} . 



This expression is obtained using Lemma [4.1| . Now it is clear that U a maps H s 
vector fields into H s vector fields since H s ~ 2 forms a multiplicative algebra, and 
since the fluctuation tensor F at t = is given by Fq which is C°°. In particular, in 
the Lagrangian frame, F is frozen, so the elliptic operator [(1 — a 2 A)(r) o n^ 1 )] o rj 
has C°° coefficients. 



Thus, the proof of Theorem 2 in Shkoller 2000 | gives a unique curve n € 
C°°(I,'Df lD ) solving (4.4). That u is only C° in time follows from the fact that 



the map n n' 1 : -> V S ^ D is only C°. That F is in C°(I, iP-^T 2 ' )) follows 
from the regularity of n. ■ 



4.2 A Corrector for the Macroscopic Velocity 
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4.2 A Corrector for the Macroscopic Velocity 



The solution to the anisotropic averaged Euler equations (|3.7| ) and (|3.8|) yields the 
pair (u,F). The macroscopic spatial velocity field u is only the zeroth-order term 
in the expansion (in e) for the velocity field u e ' S . We have computed, in equation 
(|2.6| ), the expansion of u e,d to order 0(e 2 ) as 



u e ' e {t,x) = u(t,x) +2eDef(u)(t,x) -£'(t,x,9) + 0(e 2 ). 

Since |e| is bounded by a/2, we have that 

u a > e (t, x) = u(t, x) + a Def (u)(t, x) • f (i, x, 9) + <3(a 2 ), 

so that we may add the 0(a) term to the expansion by solving for the infinitesimal 
fluctuation vector £'(t,x,9). This, however, only requires the solution of the simple 
linear advection problem (|2.5|) given by 



i'(t,x,6) b + £ uit!x) £(t,x,9) b = 0. 

Computationally, this means that we may solve for the macroscopic velocity field 
u at spatial scales larger than a and correct for the unresolved small scales to 0(a 2 ). 



4.3 Limits of Zero Viscosity 

Pcskin [1985J showed that by perturbing the Euler solution's Lagrangian particle 
trajectory rjt(x) by Brownian motions and averaging over such perturbations, the 
Navier-Stokes equations are obtained. In other words, letting Euler trajectories take 
random-walks produces the viscosity term uAu, where rjt is the flow map for the 
velocity field u. In the setting of the averaged Euler equations, the Lagrangian 
trajectory rj t (x) of a particle x corresponds to the flow of the velocity u(t, x) solving 
the anisotropic averaged Euler equations. Thus, Peskin's argument can be carried 
over in this setting to obtain the same viscous term vAu. 

We are hence motivated to define the anisotropic averaged Navier-Stokes 
equations by 

dt(l ~ a 2 C)u y + vV(l - a 2 C)u y - a 2 [Wf • Cu v + 2a 2 F : 

-4a 2 Div ((Defu") 2 ■ f) = - gradp + uAu,", u>0 (4.5) 

together with the advection for the fluctuation tensor F given by ( |3.8| ), the incom- 
pressibility constraint divn = 0, initial data n(0) = no and F(0) = Fq, and the 
no-slip conditions u = on the boundary. 

Let r\i denote the Lagrangian flow of the solution u v of the anisotropic averaged 
Navier-Stokes equations ( f4.5[ ), and let if denote the partial time derivative of the 
flow, i.e., the material velocity field. 

Theorem 4.4. For s > (n/2) + 2 and (u ,F ) € V* x C°°(T 2 ' ), there exists a 
T > 0, depending only on \\uo\\h s on not on the viscosity v, such that for each v > 

t) isinC°°([0,T),TZ^ D ), 



V(Def^ b ) 2 
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and has C°° dependence on the initial velocity field uq. Furthermore, u v t 
is in C°([0, T),V„) D C r ([0, T), V*~ r ) and depends continuously on uq. 



Vt °V t 



The proof follows the proof of Theorem 2 in Bhkoller 200C ] ; we refer the inter- 
ested reader there for the details. As a consequence of the time interval [0, T) of 
solutions u u being independent of v, we immediately obtain the following. 



Corollary 4.5. For s > (re/2) + 2, solutions u u of (4^5) converge regularly to the 



inviscid solutions u of ^3. T j as v — > 0. Furthermore, letting u v = dtr\ v o tj~ , the 
viscous Lagrangian flow rj v converges regularly in the H s topology to the inviscid 
Lagrangian flow i] = rp . 

This result states that we can generate smooth-in-time classical solutions to the 
anisotropic averaged Euler equations by obtaining a sequence of viscous solutions 
and allowing v to go to zero, and what is surprizing, this holds even in the presense 
of boundaries. Results of this type were conjectures in Marsden, Ebin and Fischerj 
1972 1 and Barenblatt and Chorin 1998a| (see also Barenblatt and Chorin [1998b|). 



In the case of the isotropic averaged Euler equations, |Foias, Holm, and Titi 
2000| have added the dissipative term z/A(l — a 2 A)u instead of using uAu, and 
this is enough to give global in time classical solutions in dimension three. It is, 



however, the term uAu that arises from either the approach of Peskin [1985 ] noted 
above, or from the constitutive theory approach of [Rivhn and Erickson 1955], 



Future Directions 

There are several interesting directions 



1. Of course numerical implementation for specific flows will be of great interest. 

2. Modeling the mean velocity profile for turbulent flows in channels and pipes. 

3. Specific flows and special solutions. 

4. Links with elliptical vortex blob methods of Zabusky and coworkers (see, e.g., 



Mclander, Zabusky and McWilliams 1988| ) would be of interest to establish; it 



is reasonable to expect that solutions of this sort would provide the anisotropic 
analog of the vortex blob solutions of Oliver and Shkoller 2000| 1 . 



5. Further investigation of the vorticity formulation and its relation with the 
coadjoint orbit structure in the semidirect product for the Hamiltonian version 
of this theory. 
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